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Measurement of the LT-asymmetry in tt° electroproduction at 
the energy of the Z\(1232) resonance 
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Abstract. The reaction p(e,e'p)ir° has been studied at Q 2 =0.2 (GeV/c) 2 in the region of W=1232 MeV. 
From measurements left and right of q, cross section asymmetries plt have been obtained in forward 
kinematics plt(Q c ^S = 20°) = (-11.68 ± 2.36 sta t ± 2.36 s „ s ) and backward kinematics plt{0 c ^ = 160°) = 
(12.18 ± 0.27 stat ±0.82 s „ s ) tt°. Multipole ratios SR{S' 1+ Mi + }/|A/i + | 2 and %i{SS + M 1+ }/\M 1+ \ 2 were de- 
termined in the framework of the MAID2003 model. The results are in agreement with older data. The 
unusally strong negative ^R{Sq + Mi+}/\Mi+ | 2 required to bring also the result of Kalleicher et al. in ac- 
cordance with the rest of the data is almost excluded. 



PACS. 13.60.Le Meson production 
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1 Introduction 

The nucleon ground- and excited state properties pres- 
ently elude a consistent description in terms of QCD as 
the basic theory of strong interaction, due to the non- 
linear, non-perturbative interaction of quarks and glu- 
ons. Over the last years, considerable efforts aimed at a 
better understanding of this complicated structure, both 
theoretically and experimentally One important issue is 
the understanding of the 'shape' of the nucleon. Despite 
its spin of 1/2 and, in consequence, the vanishing spec- 
troscopic quadrupole moment, the nucleon wave function 
might have quadrupole components which are expected to 
exhibit in the transition of the ground state to the spin 
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3/2 Z\(1232) excitation. Within constituent quark mod- 
els those components originate from tensor forces gen- 
erated by a color hyperfine interaction [l]|2]|^Pj. Larger 
quadrupole strengths are expected from models emphasiz- 
ing the particular role of pions via exchange currents [5] 
or the 'pion cloud' (i 7 S Oil 101. and also in first quenched 
Lattice QCD calculations . Dynamical approaches [111 
I12H13| enable a decomposition into the "bare" contribu- 
tions, as described in quark models, and the "dressing" by 
the pion cloud. 

The quadrupole strength is usually characterized by 
the ratios Rem — E 1+ /M 1+ and i?sM = S\ + /Mi + of 
the ttN multipoles in the Z\(1232) — > Nir decay 1 , which 
are uniquely related to the photon multipoles of electro- 
magnetic excitation |17II18| . Hence, these ratios can be 
measured in photo- and electroproduction of pions in the 
energy region of the ,4(1232) resonance. Since unwanted 
non-resonant contributions are strongly suppressed in the 
7T° channel compared to the charged pion production, most 



measurements focused on the — + pn° reaction. 

A number of studies pursued at the laboratories pro- 
viding cw electron beams yielded precise coincidence data 
based on high luminosity beams and high resolution de- 
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tectors with large angular coverage. Partially single or 
double polarization observables have been measured. The 
evolution of -Rem and Rsm with negative squared four- 
momentum transfer, Q 2 , has been investigated over a large 
range in Q 2 up to 4 (GeV/c) 2 [T$ll2tjll2T] . The extraction of 
the quadrupole ratios from the measured cross sections is 
non-trivial. For the Rsm discussed here, it is more reliable 
at lower Q 2 where M\ + dominance is more pronounced 
than at higher Q 2 and single |22U23j and double polar- 
ization results 25 ,26 are already available in addition 
to unpolarized recent measurements |20II21II27II28II29II30) . 
The low Q 2 results are almost all compatible with each 
other, yielding Rsm — —6 %, cf. fig. El The only exception 
is the result of Kalleicher et al. 27 . However, due to the 
particular kinematics it could be interpreted in line with 
the other results, if the ratio S 0+ /M 1+ ~ —10 % [THJ. #0+ 
is related to the spin 1/2 — > 1/2 transition. This amplitude 
was neglected in the analysis of • Both magnitude and 
sign of such an So+ are however unexpected from mod- 
els, e.g. MAID2003 18,, but not excluded by older mea- 
surements with large errors [3HI32] which yielded slightly 
positive values with errors of the order 10 % absolute. 

In order to investigate this issue, measurements of ir° 
electroproduction in forward and backward direction have 
been performed, which are reported in this paper. It is or- 
ganized as follows: In the next section the cross section 
formalism is briefly summarized and the method is moti- 
vated. The description of the experiment is then followed 
by a discussion of the data analysis, systematic error con- 
tributions and the results in sections El and [7| 



2 Cross section of pion electroproduction 

In one-photon-exchange approximation the fivefold differ- 
ential cross section of pion electroproduction 
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factorizes into the virtual photon flux 
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(1) 



(2) 



and the virtual photon cm cross section d 2 a v /dfi'^ n . 
Here a denotes the fine structure constant, fc 7 = (W 2 — 
m 2 )/2m p the laboratory energy of a real photon for the 
excitation of the target with mass m p to the cm energy W, 
and e = [l + (2|q| 2 /Q 2 ) tan 2 -rf]" 1 the photon polarization 
parameter. Q 2 = \q\ 2 — uj 2 is the negative squared four- 
momentum transfer, q and lu are the three-momentum 
and energy transfers, respectively, and E, E' and i? e the 
incoming and outgoing electron energy and the electron 
scattering angle in the laboratory frame. 

The unpolarized cross section for pion production with 
virtual photons is given by [171118] 



dm 



=: a v = a T + e L a L 



The partial differential cross sections, for which we use 
the short-hand notation Uj, describe the response of the 
hadronic system to the polarization of the photon field, 
characterized by the degrees of transverse (T) and longitu- 
dinal (L) polarization, e and cl = e, respectively. The 
angle <fi is the tilting angle between the electron scattering 
plane and the reaction plane. At 4> = 0° and 180° (<fi = 
90° and 270°) pions are ejected in (perpendicularly to) the 
scattering plane. 



3 Method 

The partial cross section ult is sensitive to both So+ and 
Si+. It can be determined from a fit of the (^-dependence of 
the cross section of eq. (PJ- To this end, two measurements 
left (4> = 0°) and right (4> = 180°) of the q direction are 
sufficient, which allow to form the asymmetry 



PLT(e c J) 



= Q°)-a v ((j}= 180°) 
= 0°) + <7„(<£ = 180°) 



(4) 



as a function of the 7r° center-of-mass polar angle, 9%q. 
According to eq. it is related to the partial cross sec- 
tions via 



y/2e L (e + 1)<T LT 
or + + e<JTT 



(5) 



The sensitivity to So+ and S\+ is shown by a partial wave 
decomposition of eq. ©, where only the leading multi- 
poles are retained. At the Z\(1232) resonance position the 
asymmetry 



/ ncm \ 



MrO ) -p (6) 



is obtained. Thus measurements of Plt in the forward 
(6*i ) and backward cm-hemisphere (82= n - 61) allow the 
extraction of 5i+/Mi+ and Sq+/Mi + : 



$i{S* 1+ M 1+ } 

\M 1+ \ 2 
9K{S » + M 1+ } 
\M 1+ \ 2 



/i(0i, 2 ) • [plt(0i) - Plt(0 2 )} + Ci (7) 



= MO1.2) ■ [ P lt(0i) + plt(8 2 )} + C . (8) 



+ x/2ez,(l + e)(JLT coscf) + earr cos 20. (3) 



The functions f 0(81,2) and / 1(6*1,2) denote kinematical fac- 
tors, Co and C\ contain contributions of multipoles be- 
yond the approximation. 



4 Experiment 

The p(e, e'pjir experiment was performed at the Mainz 
Microtron MAMI using a beam energy of 855 MeV 
and currents of ~ 33 \xA which were measured with high 
precision by a Forster probe in the recirculation path of 
the 3 rd microtron stage. The beam hit a liquid hydro- 
gen target. Specifically designed for this experiment, the 



D. Eisner, A. Siile et al.: Measurement of the LT- asymmetry in tt° electroproduction at the energy of the Z\(1232) 3 



Kin I Kin II Kin Ma 




100 -50 50 100 -100 -50 50 100 -100 -50 50 100 

m miss -ny (MeV/c 2 ) m miss -m„. (MeV/c 2 ) m miss -nv (MeV/c 2 ) 



Fig. 1. Typical coincidence-time (a,b,c) and missing-mass (a',b',c') spectra for Kin. I, II and Ha. Light spectra result from 
standard cuts without phase-space restrictions. Shaded time spectra (FWHM peaks: (a) 0.8 ns, (b) 2.7 ns, (c) 3.0 ns) result 
from missing-mass cuts (dashed vertical lines in the missing-mass plots) around the ty° mass; in Kin. I and la (not shown) 
the cut eleminates the 7r~-time peak and at Kin. II and Ha the prompt time peak becomes symmetric (see text). The shaded 
missing-mass spectra result similarly from the indicated cuts around the coincidence-time peak. 



1 cm cylindrical target cell with 6.25 fim Havar walls 
enabled the detection of very low-energetic protons. 
The scattered electrons were detected at a central angle 
of 0^L 6 =44.45 C and central momentum of p=408.7 MeV/c 
in Spectrometer A of the Three-Spectrometer setup of 
the Al-collaboration [35]. It consists of a QSDD mag- 
netic system and is equipped with two double planes of 
vertical drift chambers for measurement of particle tra- 
jectories in the focal plane. During the course of the mea- 
surements presented here, the standard Cherenkov detec- 
tor for 7r _ /e _ -discrimination was not available, since it 
was replaced by a focal-plane proton polarimeter |36J for 
other experiments |24U87| . In coincidence with the scat- 
tered electron, the recoil protons of the p(e, e'pjw re- 
action were detected in Spectrometer B with a similar 
focal-plane instrumentation. The smallest possible angle 
between Spectrometer B and the exit beam-pipe is 9° and 
the momentum-threshold for the proton-detection is 250 
MeV/c. Hence Q 2 =0.2 (GeV/c) 2 was the minimum possi- 
ble momentum transfer that could be reached at W=1232 
MeV. The four different kinematic settings are summa- 



Table 1. Proton kinematics. 



Kin. 


e c j n 


<M°) 


p l p ab (MeV/c) 


e l p ab (°) 


I 


160 





741.7 


33.0/32.0 


la 




180 




20.9/21.9 


II 


20 





265.02 


44.2/43.7 


Ha 




180 




9.8/10.3 



possibly caused by inefficiencies of the focal plane detec- 
tors in the proton arm, Spectrometer B was displaced by 
1° against the nominal setting for part of the measure- 
ments. 

The 7r°-data were supplemented by elastic p(e, e'p) mea- 



Table 2. Elastic scattering kinematics. 



Spec. 


A (electron) 


Spec. B (proton) 


Beam 


on 


p (MeV/c) 


9 (°) p (MeV/c) 


E (MeV) 


55.5 


612.0 


43.6 - 46.0 704.4 


855.0 


46.5 


314.6 


58.4, 59.4 251.0 


351.3 



surements (table to monitor the overall experimental 
consistency with high precision. The overdetermined kine- 
matics allows comparison of every measured variable with 
the values calculated from the other measured variables. 
Corrections of 0.5 MeV/c for the central momentum of the 
electron-spectrometer and 0.1° for the angle of the proton 
spectrometer were determined. The probable origin is a 
very slight mismatch between hardware (detector angle, 
field integral) and the track reconstruction. 

Precise measurements of electron beam current and 
dead time allowed an accurate determination of the effec- 
tive luminosity. 



5 Data analysis 



rized in tabled In order to check for false asymmetries, 



Typical coincidence-time and missing-mass spectra are 
shown in fig. ^ The overdetermined kinematics allows 



4 D. Eisner, A. Siile et al.: Measurement of the LT-asymmetry in ir° electroproduction at the energy of the Z\(1232) 



180 




1300 



nlab 

v P - 



33° 



and 



1150 1200 ^,^ 1250 

W (MeV) 

Fig. 2. Spectrometer acceptances for left 
right (0 l p ab —20° ) settings of kinematics I. W is limited by the 
electron-spectrometer (A) acceptance. The angular acceptance 
of the proton-spectrometer (B) determines the width in Bffi. 



the reconstruction of the unobserved ir° by its missing 
mass, m 1T rniss . Basic background reduction is obtained by 
coincidence-time cuts and subtraction of random coin- 
cidences via sidebands. In addition to the almost back- 
ground-free e'-p coincidence peak, the time spectrum for 
the high proton-momentum kinematics shows a smaller 
second peak at ~ -2.2 ns (Kin. I, cf. fig.^). It is caused by 
negative pions, predominantly from 7r + 7r~ reactions, the 
7r~ of which are detected in the electron spectrometer after 
a longer flight time compared to electrons. These events 
can be eliminated by the coincidence-time cut indicated 
in fig.^i. However, for Kin. II and Ha, the unwanted neg- 
ative pions can no longer be separated by coicidence time, 
due to insufficient time resolution caused by multiple scat- 
tering at the low proton-momentum. Instead, additional 
missing-mass cuts are used to suppress these events. As 
also illustrated in fig.m with a cut around m w a (a', b', c') 
the 7r~ peak vanishes in Kin. I (shaded area of fig. |TJi). Un- 
der the conditions of Kin. II/IIa the resulting coincident- 
time peak becomes symmetric after the missing-mass cut. 

Standard cuts ensure valid track reconstruction in both 
spectrometers. No target- vertex cuts were applied in order 
to avoid artificial pLT-asymmetries from the very differ- 
ent vertex-resolution along the beam-direction for the dif- 
ferent settings. Spectrometer acceptances were normalised 
with standard Monte-Carlo phase-space simulations, 
which also include the radiative corrections . 

The limited spectrometer acceptances cause different 
correlations between W, Q 2 ,e,0!£o and (f> for the settings 
left and right of q, as illustrated in fig.|2l Due to these cor- 
relations, equal binning in the variables nevertheless leads 
to unequal distributions left and right of q. Thus artificial 
PLT-asymmetries can be generated, if the mean values of 
the kinematic variables differ between left and right. This 
is obvious, e.g., for a case where Wi e ft — 1232MeV—S and 
W r ight = 1232MeV + S, since the trivial W dependence 
of the cross section produces a plt 7^ 0. 

It is extremely important to base the experimental 
asymmetries on left-right bins with the same mean values 
of the variables W, Q 2 , 6*£™ and </>. This is ensured by pro- 
jection of the numbers measured in each bin to the same 



"nominal kinematics". For this projection we made use 
of the MAID2000 parametrisation. The projection factors 
are obtained as the calculated ratios of differential cross 
sections. In order to minimise the projection error, only 
data are used within the overlap region of the two 

acceptance bands in fig. [2] Remaining uncertainties are 
included in the systematic error. 

The appropriately normalised and projected numbers 
of events left (1) and right (r) of q are determined by 



n l(4>=0°),r(4>=180°) 



E B in S ( N l,r/Pl,r)-MATD% 



(9) 



Here Ni_ r denotes the number of counts after cuts, which 
has to be divided by the relative phase-space acceptance 
Pl >r . MAID^ rr is the projection factor and L; jr repre- 
sents the relative luminosity. The asymmetry is then sim- 
ply given by 



Plt(9 c J) = 



n\ — n r 



(10) 



6 Systematic errors 

The systematic error has been estimated for Kin. I/Ia from 
the data themselves by variation of all kinematic cuts. For 
the data with low proton-momentum (Kin. II/IIa) such an 
analysis is limited by the available statistics and experi- 
mental resolution. The sliding cuts in the variables W and 
m miss resulted in non-negligible systematic errors (table 
EKjEI 3 )- The sliding cut in ra r \ liss sets a limit on remaining 
radiative effects beyond those included in the phase-space 
simulation. The spectrometer correction, which was de- 
termined by elastic measurements, has been taken into 
account both in analysis (track reconstruction) and simu- 
lation. The value given in table [3]: results from the varia- 
tion of the angle of Spectrometer B by ±0.1°. Potential er- 
ror contributions of the MAID-projection were estimated 
through relative variation of M\ + by ±5% and, simulta- 
neously, of Si+ and S 0+ by ± 50% in the full MAID2000 
calculation. The largest deviation is given in table|3Ji. Ad- 
ditional errors for the luminosity determination are not re- 
quired. The maximum variations of 2%-relative were cor- 
rected, and the remaining effect is negligible. All kine- 
matic settings were measured repeatedly to avoid time- 
dependent effects, e.g. efficiency variations. These data 
sub-sets were combined for the left- and right- kinematics. 



Table 3. Absolute systematic errors (Aplt) 
momentum kinematics 

a) W cut 0.29 % 

b) rri£ iss cut 0.23 % 

c) Spectrometer corrections 0.57 % 

d) MAID2000 projection 0.46 % 



of high proton- 
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Fig. 3. The measured plt asymmetries compared with model 
predictions from MAID2003 .18 (dotted), DMT2001 |llll2| 
(dashed), Sato/Lee |13| (dashed dotted). The full curve repre- 
sents the MAID2003 re-fit reported in this paper. The depicted 
errors represent the statistical (inner bars) and the quadratic 
sum of the statistical and systematical errors (outer bars) as 
discussed in the text. 



7 Results and discussion 

From eq. EH the asymmetries 

Plt{0 c J = 160°) = (12.18 ± 0.27 stat ± 0M sys )% 
Plt{0 c J = 20°) = (-11.68 ± 2.36 stat ± 2M sys )% 

are determined. The total systematic error is obtained by 
quadratic summation of the individual contributions in ta- 
ble [31 For the forward measurement a systematic error of 
the same size as its statistical one is assumed as a worst 
case estimate. Using these new data in conjunction with 
the previous measurement of the plti asymmetry (fifth 
structure function) of Bartsch et al. |53], we performed a 
re-fit of the MAID2003 parameters. We obtained sensitiv- 
ity to real and imaginary parts of the Si+ and So+ am- 
plitudes in the pir° channel. The results for plt and plt/ 
are depicted in fig. [3] and 0] which, for comparison, also 
shows the standard MAID2003 and the calculation within 
the dynamical models of Kamalov/Yang (DMT2001) [TT1 
E] and Sato/ Lee [13]. 



Table 4. Comparison of multipole ratios from data and cal- 
culations, as discussed in the text. 





R{5r+Af 1+ > ^ 
|M 1+ | 2 y /o > 


M{So+ M i+} toy \ 
|M 1+ |- ( /0 > 


MAID2003 re-fit 


-5.45±0.42 


2.56±2.25 


from eqs. (17151 


-4.78±0.69 


0.56±3.89 


MAID2003 


-6.65 


7.98 


Sato/Lee 


-4.74 


5.14 




-0.08 - 



-0.1 - 



120 130 140 150 160 170 18C 

CWeg) 

Fig. 4. Results for plt 1 from reference |23| with model 
predictions from MAID2003 18. (dotted), DMT2001 jllllS) 
(dashed), Sato/Lee |13| (dashed dotted). The full curve repre- 
sents the MAID 2003 re-fit. The depicted errors are only sta- 
tistical. 



From our MAID re-fit we extract the results given in 
the first row in table Q] The denoted errors are due to the 
re-fit of S\+ and So+ within the MAID2003 analysis taking 
into account the statistical and systematical errors. The 
model dependence of the extraction can be estimated from 
the truncated multipole result given in the second row in 
tabled] In the framework of this approximation we extract 
from the measured plt asymmetries the multipole ratios 
via eqs. (0and|HJ) with only leading terms in 5o+ and 
Mi+. The last two lines in table 0] contain standard model 
values without re-fit to our data. 

Figure El shows our full MAID2003 result for R SM . 
Within the errors, our extracted value is in accordance 
with measurements at the same Q 2 [2U and measurements 
at adjacent Q 2 |2Tll2^] . The negative-slope tendency of 
the CLAS-data [21] seems to be further supported by our 
Rsm value at smaller Q 2 . Provided there is no sharp Q 2 - 
dependence in Rsm, we can rule out the result of Kalle- 
icher et al. [27J at Q 2 =0.127 (GcV/c) 2 . This has been 
argued before from measurements in backward kinemat- 
ics |29| where an Sq+ contribution can not be excluded. 
In contrast, our conclusion comes from forward kinemat- 
ics as also exploited by [21] • This is illustrated in in fig. 
E] On a different scale it shows the same p^T-asymmetry 
as fig. [3 the Maid2003 re-fit and a full MAID2003 calcu- 
lation using the multipole ratios Sq+/AIi+ and S\+/M\+ 
of Kalleicher et al.. In addition, a full MAID2003 calcula- 
tion using 5o+/Mi + ~ —10%, which could reconcile the 
Kalleicher result with others |lrij . is shown. However, this 
is clearly excluded by our measurement at d c ™ — 20° with 
a 4cr significance. 

Our S Q+ /M 1+ ratio, extracted from the MAID2003 re- 
fit, is plotted in fig. 0Within the errors it agrees with 
older data at slightly larger Q 2 [3*T1I52"| . Although a little 
lower, it is also compatible with the Sato/Lee, DMT2001 
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Fig. 



Q 2 

5. Result for 



with statistical and 



systematical errors as extracted from this experiment using 
the MAID2003 re-fit (full cross), compared to measurements. 
Data where only statistical errors are given: DESY |31| (open 
square), NINA 32 (open circles), Bonn synchrotron 1 101 (open 
triangle tip up) and ELSA |57j (full triangle tip down) . Data, 
where statistical and systematical errors are given: ELSA |28| 
(full circle, to improve the presentation shiftet from Q 2 =0.201 



(GeV/c) 2 to Q 2 =0.221 (GeV/c) 2 ), MAMI 



(open dia- 



mond), CLAS ,21 (full triangles) and BATES (full 
square). The curves show model calculations MAID2003 |18| 
(solid), DMT2001 [12] (dashed) and Sato/Lee [E] (dashed dot- 
ted). 



and standard MAID2003 parametrisations. In view of the 
quite large experimental errors it is not yet clear whether 
this ratio differs from zero. But we can not support a large 
negative 5 + /Mi+ ratio. 

A sensitive access to the ratio S0+/S1+ is provided by 
a precise measurement of the zero-crossing of plt or <tlt- 
By now it is possible to extract this ratio from available 
data close to the zero crossing at Q 2 —Q.121 (GeV/c) 2 |3U] . 
with the result compatible to the MAID2003 parametri- 
sation. Other existing data [23 cover the full range of 8",] 1 
at Q 2 =0.4-1.8 (GeV/c) 2 . However, at higher Q 2 the S Q+ 
extraction seems to be affected more strongly by higher 
partial waves than expected in the paper by Joo et al. 
|21j . This might be resolved by very recent polarisation 
data [TTTFT^ . 

In future experiments at MAMI-C a more accurate de- 
termination of Sq+ at low Q 2 is feasible, using the Three- 
Spectrometer setup of the Al-collaboration complemented 
by the KaoS-spectrometer [4*3"] . 



0.4 
0.3 
0.2 
0.1 

Q. 

-0.1 

-0.2 
-0.3 
-0.4 





1 


1 1 1 1 1 1 




1 1 1 z 


r 






/' 










/ 


\ 






/ 




\ J 






/' 




\ ; 






/ ^ 




^x-N - 






















;\ 




/ / 






r n 




/' / 








\ 


/' / 








\\ 










/ / 








\\ 


// 








\';. 


--. ..<?'■'"' 








1 







1 , ] , = 



20 



40 



100 
er'(deg) 



80 

.cm 



120 140 160 180 



Fig. 6. The plot shows, at smaller scale, the MAID2003 
re-fit (solid) in comparison to the full MAID2003 calcula- 
tion with modified ^{51+} and Q ! {S'o+} for two situations: 
5i + /Af 1+ = -12.5%, S 0+ /M 1+ = 0.0% (dashed dottet), i.e. 
the result of E3, and Si+/M 1+ = -10%, S +/M 1+ = -14% 
(dottet), which would bring [27) in accordance with other mea- 
surements |16) . 
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Fig. 7. Result for Sft{5o + Mi + }/|Mi+ 1 2 with statistical and 
systematical errors as extracted from this experiment using 
the MAID2003 re-fit (full cross); compared to measurements 
from DESY |31| (open square), NINA |32[ (open circles), where 
only statistical errors are indicated. The curves show model 
calculations MAID2003 [THj (solid), DMT2001 Q3 (dashed) 
and Sato/Lee ^5] (dashed dotted). 



8 Summary 

We have measured the plt asymmetry in forward 
(6)=™=20 o ) and backward ((9^=160°) kinematics of ir° 
electroproduction off the proton at Q 2 =0.2 (GeV/c) 2 
around W=1232 MeV. The measurement of the two kine- 
matic settings allows the extraction of Si+ and So+ hi a 
very transparent way within a simple s- and p-wave ap- 



proximation or, alternatively, using the full MAID2003 
parametrisation without any truncation. Our results for 
S1+/M1+ and 1S0+/M1+ are in agreement with existing 
measurements and calculations. Our result removes a re- 
maining possibility to reconcile the datum of Kalleicher et 
al. for the ratio S1+/M1+ with other measurements 
through a large negative Sq + /M\ + . 

We thank T.-S.H. Lee and T.Sato for providing their cal- 
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